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ABSTRACT 

Using  numerical  models  of  force-free  magnetic  fields,  we  have 
examined  how  the  shearing  of  footpoints  in  arcade  geometries  leads  to  an 
inflation  of  the  coronal  magnetic  field.  For  each  of  the  shear  profiles 
considered,  all  of  the  field  lines  become  elevated  compared  with  the  potential 
field.  This  includes  cases  where  the  shear  is  concentrated  well  away  from  the 
arcade  axis,  such  that  Bz,  the  component  of  field  parallel  to  the  axis,  increases 
outward  to  produce  an  inward  Bz2/8n  magnetic  pressure  gradient  force.  These 
results  contrast  with  an  earlier  claim,  shown  to  be  incorrect,  that  field  lines 
car.  sometimes  become  depressed  as  a  result  of  shear.  We  conjecture  that  an 


inflation  of  the  entire  field  will  always  result  from  the  shearing  of  simple 
arcade  configurations.  These  results  have  implications  for  prominence 
formation,  the  interplanetary  magnetic  flux,  and  possibly  also  coronal  holes. 
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I.  INTRODUCTION 

A  majority  of  observed  solar  phenomena  are  either  caused  by  or 
significantly  influenced  by  the  magnetic  fields  which  permeate  the  solar 
atmosphere.  Magnetic  fields  give  structure  to  the  atmosphere;  they  direct  the 
flow  of  mass  and  many  forms  of  energy;  they  may  be  an  important  source  of 
coronal  heating;  and  they  power  many  of  the  most  spectacular  solar  events, 
most  notably  solar  flares.  It  is  not  surprising,  therefore,  that  considerable 
effort  has  been  devoted  to  studying  the  properties  of  solar  magnetic  fields. 

Deep  within  the  photosphere  and  below,  the  magnetic  field  is  dominated 
by  the  highly  massive  and  reasonably  highly  conducting  gas  in  which  it  is 
embedded.  Any  motions  of  the  gas  will  cause  the  magnetic  field  to  be  dragged 
about  the  solar  surface.  In  the  overlying  corona,  on  the  other  hand,  the  gas 
density  and  pressure  are  sharply  reduced  and  magnetic  forces  tend  to  greatly 
exceed  those  associated  with  the  material.  In  this  situation  the  Lorentz  force 
must  very  nearly  vanish  (for  equilibrium),  and  it  is  therefore  useful  to  study 
the  properties  of  so-called  force-free  magnetic  fields,  for  which 

(VxB)xB  =  0  .  q  X) 

One  important  property  of  these  fields  leir  energy  content.  As 
photospheric  and  subphotospheric  motions  displace  magnetic  footpoints, 
electric  currents  are  induced  in  the  corona.  The  energy  associated  with  these 
currents  is  known  as  the  "free  magnetic  energy"  and  is  energy  which  can 
libera  leu  from  the  field  in  the  form  of  a  flare,  for  example.  Klimchuk, 
Sturrock,  and  Yang  11988)  have  recently  shown  that  systematic  shearing 
motions  in  arcade  configurations  will  produce  energy  increases  that  can  be 
simply  and  prediclanly  described. 
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Another  important  property  of  force  free  magnetic  fields  is  their 
structure,  and  in  particular,  the  manner  in  which  this  structure  changes  as 
the  footpoint  positions  are  changed.  More  than  20  years  ago,  Sturrock  and 
Woodbuiy  (1967)  demonstrated  that  magnetic  field  configurations  can  become 
inflated  as  the  footpoints  are  sheared;  that  is,  the  field  lines  rise  and  spread 
out  horizontally.  A  number  of  later  calculations  revealed  this  to  be  a  common 
behavior  (e.g.,  Barnes  and  Sturrock  1972;  Low  1977;  Jockers  1977;  Birn, 
Goldstein,  and  Schindler  1978;  Priest  and  Milne  1980;  Yang,  Sturrock,  and 
Antiochos  1986;  Zwingmann  1987;  Mikic,  Barnes,  and  Schnack  1988; 
Klimchuk,  Sturrock,  and  Yang  1988;  Biskamp  and  Welter  1989;  Klimchuk  and 
Sturrock  1989;  Finn  and  Chen  1989;  Amari  et  al.  1989).  This  is  not  an 
unexpected  result,  since  shearing  motions  increase  the  strength  of  the  field, 
and  the  associated  increase  in  magnetic  pressure,  which  is  equivalent  to  an 
increase  in  thermal  gas  pressure  (see  below),  should  cause  an  overall 
expansion  of  the  configuration. 

We  are  aware  of  at  least  one  calculation,  however,  where  an  opposite 
behavior  seems  to  occur.  In  one  of  two  cases  studied  by  Low  and  Nakagawa 
('1975),  it  was  reported  that  the  field  becomes  depressed  (i.e.,  the  field  lines 
descend )  as  the  shear  of  the  footpoints  increases.  This  is  an  intriguing  and 
potentially  very  important  result  that  has  received  considerable  attention  in 
the  review  literature  (e.g.,  Birn  and  Schindler  1981;  Priest  1982a;  Low  1982). 

There  are  a  variety  of  reasons  for  wanting  to  study  how  magnetic 
configurations  change  in  response  to  footpoint  shearing  motions.  As 
configurations  become  inflated,  for  example,  closed  field  lines  are  brought  to 
greater  heights  in  the  atmosphere,  where  they  are  more  likely  to  be  opened  up 
to  become  part  o t  the  outflowing  solar  wind  (e.g.,  Pncuman  and  Kopp  1971). 


i  m- 


pjenees  tor  i  tie  'ntcrplaretnry  magnetic  fieki,  since  a.~>  mere 
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field  lines  open,  the  quantity  of  interplanetary  magnetic  flux  increases  (see,  for 
instance,  Wang  and  Sheeley  1988).  It  is  possible  that  this  process  could  lead  to 
the  formation  or  growth  of  coronal  holes.  In  addition,  the  lateral  expansion  of 
inflating  configurations  will  tend  to  promote  flux  linkages  between  different 
parts  of  the  Sun,  since  oppositely  directed  field  lines  will  be  forced  together  and 
hence  will  be  more  likely  to  reconnect. 

The  depression  of  magnetic  field  configurations  may  ultimately  lead  to 
prominence  formation,  since  prominences  cannot  easily  form  unless  there  is 
a  dip  in  the  field  (i.e.,  a  region  that  is  concave  upward),  or  at  the  very  least  a 
region  that  is  exceptionally  flat  (see  the  discussions  in  Priest,  Hood,  and  Anzer 
1989  and  Amari  et  al.  1989).  Even  then,  the  hydrodynamics  of  prominence 
formation  is  not  well  understood  (e.g.,  Poland,  Mariska,  and  Klimchuk  1986). 
It  is  therefore  important  to  fully  understand  the  result  of  Low  and  Nakagawa 
(1975)  and  to  determine  what  types  of  situations  will,  in  general,  lead  to  shear- 
induced  depression  rather  than  shear-induced  inflation  of  force-free  magnetic 
field  configurations.  We  report  here  the  details  of  our  investigation  into  this 
problem  (Klimchuk  1989). 

II.  FOOTPOINT  SHEAR  AND  MAGNETIC  PRESSURE 

In  many  situations,  the  coronal  magnetic  field  has  an  arcade-like 
geometry  in  which  the  field  is  approximately  uniform  in  one  of  the  horizontal 
directions.  Active  regions,  for  example,  tend  to  have  two  dominant  regions  of 
opposite  polarity  separated  by  a  long,  though  not  always  straight,  neutral  line 
<or,  more  properly,  a  "weak  field  corridor"  where  the  flux  is  very  small; 
Klimchuk  198 7).  Similarly,  much  the  quiet  Sun  is  divided  into  large, 
predominantly  unipolar  regions  of  opposite  polarity  that  are  stretched  out  into 
elongated  patterns  by  a  combination  of  supergranular  diffusion,  differential 
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rotation,  and  meridional  flow  (see,  for  instance,  Sheeley,  Nash,  and  Wang 
1987).  In  both  situations,  the  overlying  corona  is  characterized  by  nested  loop 
systems  referred  to  as  arcades. 

Often,  the  length  of  the  arcade  is  substantially  greater  than  its  width, 
and  it  is  therefore  common  to  model  magnetic  field  configurations  that  are 
assumed  to  have  translational  symmetry  (i.e.,  to  be  infinitely  long).  Taking 
the  y-direction  to  be  the  vertical  direction,  with  y  =  0  corresponding  to  the 
photosphere,  and  the  z-direction  to  be  the  direction  of  invariance,  we  can 
express  the  magnetic  field  in  the  form 

Id  A  dA  \ 

B(x,y)  =  — , - -,  B  I , 

\dy  dx  I  (21) 


which  follows  directly  from  the  divergence-free  nature  of  B.  Since 


B  •  V  A  =  0  , 


(2.2) 


we  see  that  the  scalar  function  A(x,y),  which  is  simply  the  z-component  of  the 
magnetic  vector  potential,  is  constant  along  field  lines.  Contours  of  constant  A 
in  the  x-y  plane  are  simply  the  projections  of  field  lines  onto  that  plane,  and 
correspond  to  an  end-on  view  of  the  arcade. 

If  we  include  the  effects  of  gas  pressure,  the  force  balance  equation 
[equation  (l.l)j  becomes 


— (VxB)xB  -  VP  =  (). 
An 


(2.3) 
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Upon  substituting  for  B  using  equation  (2.1),  the  z-component  of  force  balance 
leads  to  the  result  that 


Bz=B7(A), 


(2.4) 


so  that  Bz  is  also  constant  along  field  lines  (e.g.,  Priest  1982a).  The  y- 
component  of  force  balance  leads  to 


— V^A  + 
4tc 


(2.5) 


which  is  the  well-known  Grad-Shafranov  equation.  It  shows  explicitly  that  the 
magnetic  pressure  associated  with  the  z-component  of  the  magnetic  field  is 
equivalent  to  gas  pressure  in  the  way  that  it  affects  the  equilibrium 
distribution  of  field  lines,  A(x,y).  Note  that  if  the  gravitational  scale  height  is 
large  compared  with  the  geometric  height,  so  that  P  does  not  depend  explicitly 
on  height,  then  it  too  is  constant  along  field  lines.  Having  established  that  Bz 
and  P  are  analogous  in  this  sense,  we  henceforth  concern  ourselves  only  with 
force-free  magnetic  fields,  for  which  P  =  0.  It  is  important  to  point  out, 
however,  that  the  choice  BZ(A)  is  not  as  flexible  as  the  choice  P(A)  from  a 
physical  point  of  view  (e.g.,  Klimchuk  and  Sturrock  1989). 

In  an  unsheared  magnetic  arcade,  which  is  a  potential  field,  all  field 
lines  lie  in  planes  of  constant  z  and  Bz  =  0.  By  shearing  the  field  (that  is,  by 
displacing  the  footpoints  of  one  or  more  field  lines  in  the  z-direction  so  that  the 
two  footpoints  of  a  given  field  line  are  separated  in  z),  we  introduce  a  finite  Bz 
component.  We  expect  the  additional  magnetic  pressure  associated  with  this 
component  to  produce  an  overall  inflation  of  the  field;  although,  in  principle, 
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some  parts  of  the  field  could  become  depressed.  This  expected  inflation  has 
been  found  to  occur  in  a  large  number  of  specific  arcade  models,  as  cited  in  the 
Introduction  [see  also  Wolfson  and  Gould  (1985)  and  Hundhausen, 
Hundhausen,  and  Zweibel  (1981)  for  a  discussion  of  inflated  fields  in  a 
spherical  geometry!. 

Suppose  that  the  distribution  of  shear  is  such  that  Bz  increases  outward 
from  the  arcade  axis  and  reaches  a  maximum  value  a  finite  distance  away 
from  the  axis.  Then  there  will  be  both  inward  and  outward  Bz2  pressure 
gradient  forces— inward  in  the  inner  part  of  the  arcade,  and  outward  in  the 
outer  part.  One  might  expect  those  field  lines  whmn  lie  underneath  the  field 
line  of  maximum  Bz  to  become  depressed,  and  those  field  lines  which  lie  above 
it  to  become  elevated.  This  would  seem  to  be  the  situation  in  which  depressed 
field  is  most  likely  to  occur,  and  is  the  situation  we  wish  to  explore. 

One  way  to  calculate  force-free  magnetic  fields  is  to  specify  the 
functional  dependence  of  Bz  upon  A  explicitly,  which  we  refer  to  as  the 
"generating  function"  method.  In  principle,  it  is  then  straightforward  to 
study  cases  where  Bz  is  maximized  away  from  the  arcade  axis.  Analytic 
solutions  are  often  difficult  to  obtain,  however,  and  only  very  recently  have 
numerical  solutions  with  this  type  of  Bz  dependence  been  calculated  (Amari  et 
a l.  1989).  Moreover,  solutions  obtained  by  the  generating  function  method  can 
sometimes  have  peculiar,  unphysical  properties,  as  discussed  in  Section  IV 
and  in  Klimchuk  and  Sturrock  (1989).  It  generally  agreed  that  it  is  physically 
more  reasonable  to  define  the  magnetic  field  in  terms  of  footpoint  boundary 
conditions,  i.e.,  the  shear,  than  it  is  to  place  global  constraints  on  Bz  (e.g., 
Sturrock  and  Woodbury  1967;  Jockers  1977;  Low  1982;  Priest  1982a, b;  Aly  1985; 
Amari  et  al.  1989;  Finn  and  Chen  1989;  Klimchuk  and  Sturrock  1989).  This  is 
the  approach  that  we  take. 
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The  relationship  between  Bz  and  shear  is  not  always  a  simple  one.  For 
instance,  it  is  not  generally  true  that  the  field  lines  with  the  greatest  footpoint 
separation  (in  z)  are  also  the  field  lines  with  the  greatest  Bz  To  see  this,  let  s 
be  the  spatial  coordinate  along  a  contour  of  A.  We  then  have  the  relationship 


dz  _  ds 

b7b. 


(2.6) 


where  =  I  VA  j  is  the  magnitude  of  the  projection  of  the  magnetic  field 
vector  onto  the  x-y  plane.  We  can  integrate  along  the  contour  to  get  the  z- 
separation  between  the  two  footpoints: 


Az  = 


B, 
B  . 


ds  . 


(2.7) 


Noting  that  Bz  is  constant  on  A,  we  have 


B*  = 


L  '  Az  , 


(2.8) 


where  L  is  the  length  of  the  contour.  Thus,  Bz  is  related  not  only  to  the 
footpoint  separation,  but  also  to  the  mean  strength  and  length  of  the  projected 
field  line.  For  a  shear  profile  in  which  Az  increases  away  from  the  neutral 

line,  it  does  not  necessarily  follow  that  Bz  will  also  increase,  since  both  /* 
and  L  ]  decrease.  How  Bz  varies  in  space  depends  on  the  details  of  the 
particular  model. 
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III.  MODEL  CALCULATIONS 

Equation  (2.1)  is  one  mathematical  representation  of  magnetic  fields 
having  translational  symmetry.  Another  makes  use  of  Clebsch  variables  a 
and  (3: 


B  =  Va  x  Vp  . 


(3.1) 


It  is  clear  that  a  and  P  are  constant  along  field  lines.  Furthermore,  it  can  be 
easily  shown  that  a  is  equivalent  to  A  in  equation  (2.1).  It  is  convenient  to 
adopt  the  form 


a  =  a(x,y),  p  =  z  -  y(x,y) ,  (3  2) 

because  then  y  measures  the  displacement  of  any  field  line  in  the  z-direction 
(Sturrock  and  Woodbury  1967).  By  specifying  a  and  y  as  functions  of  x  on  the 
plane  v  =  0  we  are  able  to  specify  both  the  normal  value  of  the  photospheric 
magnetic  field  and  the  field  line  connectivity  (i.e.,  the  relative  positions  of  the 
magnetic  footpoints)  . 

We  solve  equation  (1.1)  subject  to  these  boundary  conditions  using  a 
numerical  relaxation  technique  known  as  the  "magneto-firictional  method" 
(Yang,  Sturrock,  and  Antiochos  1987;  Klimchuk,  Sturrock,  and  Yang  1988; 
Klimchuk  and  Sturrock  1989).  Basically,  an  initial  guess  at  the  coronal  field 
relaxes  to  a  force-free  state  through  motions  that  are  proportional  to  the  local 
Lorentz  force.  The  constant  of  proportionality  can  be  interpreted  as  a  frictional 
force,  which  is  the  origin  of  the  method’s  name.  Interested  readers  are 
referred  to  the  above  references  for  more  details  (see  also  Chodura  and 
Schluter  1981  and  Craig  and  Sncyd  1986). 
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For  photospheric  boundary  conditions  on  a  we  use  the  normal  field 
distribution  ad'  jcated  by  Klimchuk,  Sturrock,  and  Yang  (1988)  as  being 
appropriate  for  solar  active  regions: 


Bv(x,0)  = 


sin 


71 

(|x|  -  c)  2 


sin 


K 


(1x1  -  c  -  d)  2 


k  0 


I  x 1  <  c 
c  <  I  x  I  <  c+1 

c+1  <  Ixl  < c+d+1  ^  -  (3.3) 

e+d+1 <  I x I  < c+d+2 

I x I  >  c+d+2 


so  that 


'  ,  4 

d  +  — 

7C 

2 

—  cos 

K 


Ttl  2 
(|x|  -  c)  2J  +  d  +  - 


X  I  <  c 


—  c  <  I  X  I  <  c+1 


a(x,0)  = 


|x| 


x  +  c-'  d  +  l  +  — 


2 

—  cos 

7C 


71 


7t 


( I x |  -  c  -  d)  2 


2 

+  — 

K 


c+1 <  1 x I  <  c+d+1 


c+d+1 <  I x I  <  c+d+2 


i x I  >  c+d+2 


(3.4) 

Here,  c  is  the  half-width  of  the  "weak  field  corridor"  found  to  separate  the 
opposite  polarity  parts  of  most  active  regions  (Klimchuk  1987),  and  d  +  2  is  the 
width  of  the  individual  strong  field  regions,  or  "plages."  We  choose  c  =  1  and  d 
=  -1 ,  in  which  case  the  units  of  our  spatial  coordinates  correspond  to  roughly 
1(F  km  on  the  Sun.  The  normal  field  distribution  By  is  shown  in  Figure  1. 


12 


To  guarantee  that  Bz  has  its  maximum  a  finite  distance  away  from  the 
arcade  axis  we  adopt  the  following  localized  shear  profile: 

y(x,0)  =  Zma<exp[-(x  -  x0)2  ,  (3>5) 

where  xq  =  4.  This  shear  profile  is  also  shown  in  Figure  1.  Bz  is  greatest  for 
the  a  =  0.62  field  line,  which  is  displaced  inward  somewhat  from  the  field  line 
of  maximum  shear  (a  =  0.5),  as  expected  on  the  basis  of  equation  (2.8). 

The  numerical  computations  are  performed  on  a  finite-difference  grid 
covering  the  domain  0  <  x  <  60,  0  <  y  <  60  (only  half  the  field  need  be  computed 
explicitly  due  to  symmetry  about  the  x  =  0  plane).  A  total  of  289  x  289  grid 
points  are  placed  nonuniformly  throughout  this  domain  to  provide  a  resolution 
of  0.11  near  the  origin  and  0.51  at  large  distances.  "Superconductor"  boundary 
conditions  (a  =  0)  at  the  outer  boundaries  x  =  60  and  y  =  60  prevent  flux  from 
leaving  the  box,  but  the  boundaries  are  far  enough  removed  that  they  have 
negligible  effect  on  the  field  near  the  source. 

Contours  of  constant  a  are  shown  in  Figure  2a  for  the  potential  field 
case  Zmax  =  0  and  in  Figure  2b  for  the  sheared  field  case  Zmax  =  6.  It  is  clear 
from  a  comparison  of  the  plots  and  from  an  examination  of  the  raw  model  data 
that  all  of  the  field  lines  are  elevated  in  the  sheared  case.  Thus,  our  initial 
expectations  that  the  inner  field  lines  might  become  depressed  were  incorrect. 

We  have  examined  several  other  shear  functions  to  determine  if  they 
also  produce  field  lines  which  are  everywhere  elevated.  These  shear  functions 
are  plotted  in  Figure  3  (note  that  y(x,0)  is  specified  only  in  the  region  of  nonzero 
flux;.  The  function  labelled  2  is  from  Klimchuk,  Sturrock,  and  Yang  (1988), 
and  the  function  labelled  4  is  repeated  from  Figure  1  for  comparison.  All  of 
the  shear  functions  have  an  amplitude  of  6,  which  is  comparable  to  the  half- 
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width  of  the  arcade  at  it?  t'ase  and  corresponds  to  a  footpoint  displacement  of  = 
6xl04  km  in  a  solar  active  region.  Force-free  field  models  were  calculated  for 
each  of  these  shear  functions  using  the  same  normal  field  boundary  condition 
as  before. 

In  each  case  all  of  the  field  lines  are  elevated  compared  with  the 
potential  field.  This  is  not  surprising  in  models  3,  2,  and  3,  where  Bz 
decreases  monotonically  outward  from  the  arcade  axis,  but  in  model  5,  Bz  is 
greatest  along  the  a  =  0.67  field  line.  The  degree  to  which  the  field  becomes 
inflated  is  indicated  in  Table  1,  where  we  give  the  change  in  height  along  the  y- 
axis  (x  =  0)  of  three  representative  field  lines:  an  inner  (a  =  0.8),  a  middle  (a  = 
0.5),  and  an  outer  (a  =  0.2)  field  line.  As  expected,  the  absolute  height  change 
increases  as  one  moves  outward  from  the  arcade  axis.  We  see  that  the 
magnitude  of  the  change  can  be  quite  substantial.  In  model  1,  for  example, 
the  outer  field  line  rises  by  an  amount  comparable  to  the  width  of  the  arcade 
itself.  Still  higher  field  lines  rise  by  even  more. 

Also  given  in  Table  1,  in  parentheses,  is  the  percentage  change  in 
height,  Ay  /  y.  We  see  that  the  heights  of  the  inner  and  middle  field  lines  have 
more  than  doubled  in  the  first  two  models.  The  manner  in  which  the 
percentage  height  change  varies  with  y  is  somewhat  different  in  the  different 
models.  In  model  1,  where  the  shear  profile  is  flat,  Bz  decreases 
systematically  outward  according  to  equation  (2.8),  and  so  too  does  the 
fractional  change  in  height.  In  models  4  and  5,  on  the  other  hand,  Bz  peaks 
midway  between  the  inner  and  middle  field  lines,  and  we  find  that  the 
fractional  change  in  height  of  the  middle  field  line  is  much  greater  than  that 
of  the  inner  field  line.  Thus,  the  field  lines  tend  to  spread  apart  preferentially 
in  the  region  of  greatest  B7?  pressure,  as  one  would  expect. 
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The  models  are  ordered  in  Table  1  according  to  the  energy  of  the  field,  U, 
given  in  the  second  to  last  column  as  the  ratio  U/Uo,  where  Uo  is  the  energy  of 
the  corresponding  potential  field  (the  same  for  all  models).  We  find,  not 
surprisingly,  that  the  fields  with  the  greatest  energy  are  also  the  fields  with 
the  greatest  inflation.  In  the  last  column  of  Table  1  we  give  the  ratio  U/Uo  That 
is  predicted  by  a  simple  formula  derived  in  Klimchuk,  Sturrock,  and  Yang 
(1988)  and  thought  to  apply  to  a  wide  range  of  arcade  configurations.  The 
agreement  between  the  predicted  and  actual  values  appears  quite  satisfactory. 

Table  2  indicates  how  the  inflation  of  a  field  varies  with  shear  amplitude 
when  the  functional  form  of  the  shear  is  held  fixed.  All  of  the  models  were 
computed  using  the  shear  function 


T<x,0)  =  Zmax 


si 


(3.6) 


and  differ  only  in  the  value  of  Zmax  (from  Klimchuk,  Sturrock,  and  Yang 
1988).  The  case  Zmax  =  6  is  identical  to  model  2  in  Figure  3  and  Table  1.  As 
expected,  the  field  becomes  progressively  more  inflated  as  the  amplitude  of  the 
shear  is  increased.  The  rate  of  inflation  is  faster  than  the  rate  of  energy 
increase  in  the  sense  that  the  fractional  change  in  height  of  the  field  lines  is 
roughly  twice  the  fractional  change  in  energy  of  the  field. 

IV.  MODEL  OF  LOW  AND  NAKAGAWA 
Our  inability  to  generate  a  model  for  which  any  of  the  field  lines  descend 
has  prompted  us  to  reexamine  the  model  of  Low  and  Nakagawa  (1975),  which 
is  thought  to  exhibit  this  behavior.  The  boundary  condition  on  a  in  that  model 
is 
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a(x,0)  =  cos  x  . 


(4.1) 


Low  and  Nakagawa  take  the  generating  function  approach  and  specify 


Bz(a)  =  -  e 


a 


1  +  a 


(4.2) 


where  e  is  a  parameter  analogous  to  Zmax.  Since  this  particular  form  of  the 
generating  function  is  not  amenable  to  analytic  solution,  the  authors  obtain  a 
numerical  solution  over  the  domain  -jt/2  <  x  <  k/2,  0  <  y  <  10,  where  it  is 
assumed  that  a  =  0  on  the  outer  boundary. 

At  this  point  we  note  an  important  property  of  equation  (4.2):  |  Bz  |  is  an 
increasing  function  of  a  and  therefore  a  decreasing  function  of  I  x  I  throughout 
the  domain  of  the  solution.  The  pressure  gradient  force  associated  with  Bz  is 
everywhere  outward.  For  this  reason  the  apparent  result  of  Low  and 
Nakagawa  that  the  field  becomes  depressed  is  suspicious. 

In  order  to  reconstruct  the  Low  and  Nakagawa  field  using  the  magneto- 
frictional  method  we  must  know  the  z-separation  of  each  pair  of  footpoints  (the 
shear  function).  This  can  be  obtained  from  the  generating  function  solution  by 
integrating  along  field  lines  according  to  equation  (2.7),  as  was  done  by  Low 
and  Nakagawa  for  five  representative  field  lines.  The  resulting  values  for  the 
case  e  =  0.2  are  plotted  in  Figure  4  (taken  from  Figure  3  in  Low  and 
Nakagawa).  In  addition,  we  know  that  the  outermost  field  line  intersecting 
the  photosphere  at  x  =  ±  n/2  has  no  shear,  because  Bz  =  0  for  this  field  line. 
Finally,  while  we  have  no  detailed  information  on  the  footpoint  separation 
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near  the  arcade  axis  (this  was  not  calculated  by  Low  and  Nakagawa),  we  do 
know  that  the  Bz  pressure  will  be  smallest  for  a  small  separation.  In  order  to 
maximize  the  tendency  for  the  field  to  become  depressed,  we  therefore  take  the 
footpoint  separation  to  vanish  at  x  =  0.  A  least-squares  fit  of  the  resulting  7 
data  points  leads  to  the  following  polynomial  expression  for  the  shear  function: 

y(x,0)  =  c4x  +  c3x  +  c2x  +C]X,  ^  ^ 

where  C4  =  -1.154,  C3  =  2.829,  C2  =  -2.543,  c2  =  1.487.  This  shear  function  is 
plotted  in  Figure  4. 

We  have  computed  a  force-free  field  by  the  magneto-frictional  method 
using  equations  (4.1)  and  (4.3)  as  photospheric  boundary  conditions  and  using 
the  same  computational  domain  and  outer  boundary  conditions  as  were  used 
by  Low  and  Nakagawa.  The  field  is  shown  in  Figure  5  together  with  the 
corresponding  potential  field  [y(x,0)  =  0J.  The  sheared  field  is  clearly  inflated 
compared  with  the  potential  field,  as  expected  on  the  basis  of  the  above 
discussion  and  in  direct  contradiction  with  the  generating  function  result  of 
Low  and  Nakagawa.  We  conclude  that  the  earlier  result  is  either  incorrect  or 
else  has  a  topology  that  is  more  complicated  than  originally  thought,  a 
conclusion  that  is  shared  by  Low  (private  communication). 

It  is  a  well-known  property  of  generating  function  solutions  that  they 
are  not  always  unique  (e.g.,  Birn,  Goldstein,  and  Schindler  1978;  Low  1982; 
Heyvaerts  et  al.  1982).  Depending  on  the  form  of  the  generating  function,  there 
may  be  two  or  more  solutions  which  satisfy  the  same  boundary  condition  on 
cdx,0)  and  have  the  same  Bz(a]  dependence,  but  which  have  a  different  field 
line  connectivity  (i.e. ,  footpoint  shear).  This  is  true  of  the  generating  function 
used  by  Low  and  Nakagawa,  for  example.  It  is  also  possible  to  have  multiple 
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solutions  which  satisfy  the  same  boundary  condition  on  a(x,0)  and  have  the 
same  field  line  connectivity,  but  which  have  a  different  Bz(a)  dependence 
(Klimchuk  and  Sturrock  1989;  see  also  Priest  and  Milne  1980).  Low  and 
Nakagawa  attempted  to  find  the  solution  having  the  simplest  topology-the  one 
for  which  all  of  the  field  lines  are  connected  to  the  photosphere.  It  may  be, 
however,  that  they  instead  found  a  solution  with  disconnected  or  "floating" 
flux  (Low,  private  communication).  If  the  quantity  of  floating  flux  increases 
with  e,  as  it  does  in  the  generating  function  sequence  of  Low  (1977),  the 
magnetic  pressure  associated  with  the  flux  might  force  the  underlying,  arcade 
field  downward,  in  agreement  with  the  Low  and  Nakagawa  result. 

V.  DISCUSSION 

The  results  of  the  previous  two  sections  lead  us  to  propose  the  following 
conjecture:  In  fields  having  a  simple  arcade  geometry  (without  floating  flux), 
any  form  of  footpoint  shear,  even  one  that  is  localized  away  from  the  arcade 
axis,  will  cause  all  of  the  field  lines  to  rise  and  the  entire  configuration  to 
inflate.  Moreover,  by  equation  (2.5),  any  enhancement  to  thermal  gas  pressure 
must  have  the  same  effect  on  magnetostatic  fields  for  which  the  pressure  does 
not  depend  explicitly  on  height.  Although  we  have  been  unable  to  prove  this  as 
a  general  theorem,  we  are  able  to  provide  two  plausibility  arguments. 

First,  Aly  (1985)  has  shown  that  as  the  footpoint  shearing  displacements 
of  an  arcade  field  become  arbitrarily  large,  the  field  approaches  the  "open 
field”  configuration  as  its  limiting  state  (although  see  Finn  and  Chen  1989). 
For  any  shear  function  Zmaxy(x)  that  is  everywhere  nonzero,  all  of  the  field 
lines  extend  to  infinity  as  Zmax  — >  •».  In  other  words,  the  field  becomes 
infinitely  inflated.  It  is  reasonable  to  expect  that  a  sheared  arcade  will 
monotonically  approach  this  limiting  state,  rather  than  first  deflating  for 
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some  range  of  shear  and  then  reversing  direction  and  inflating  for  larger 
values  of  shear.  We  see  no  reason  for  why  the  interplay  between  magnetic 
tension  and  pressure  gradient  forces  should  behave  differently  for  small  and 
large  values  of  shear. 

As  a  second  plausibility  argument  we  invoke  a  simple  hydrostatic 
analogy.  Consider  a  semi-infinite  vertical  tube  filled  with  an  isothermal, 
gravitationally  s  "ratified  gas.  The  tube  is  capped  at  the  bottom,  but  not  at  the 
top.  Suppose  that  the  temperature  of  the  particles  in  some  layer  part  way  up 
the  tube  is  increased,  as  indicated  in  Figure  6.  If  the  gas  responds 
adiabatically  to  the  temperature  increase,  a  new  equilibrium  will  be  reached 
in  which  the  hotter  gas  and  the  cooler  gas  above  are  elevated,  but  the  cooler  gas 
below  is  unchanged  (for  uniform  gravity).  This  result  comes 
straightforwardly  from  the  fact  that  the  pressure  at  any  point  in  the  gas  must 
be  equal  to  the  weight  of  the  overlying  material: 


The  weight  of  the  material  is  of  course  unaffected  by  the  temperature  increase. 
Initially  there  will  be  a  pressure  increase  in  the  heated  layer,  and  no  doubt  the 
transient  response  will  be  to  force  the  underlying  material  downward,  but 
eventually  the  atmosphere  will  adjust  such  that  none  of  the  material  is 
depressed  in  the  final  static  equilibrium.  In  fact,  if  gravity  decreases  with 
height,  as  it  does  in  the  solar  atmosphere,  all  of  the  material  will  become 
elevated  compared  with  its  original  position,  including  the  material  that  lies 
below  the  heated  layer  (as  pointed  out  by  P.  A.  Sturrock). 
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This  result  is  analogous  to  our  result  for  magnetic  fields  in  which  all  of 
the  field  lines  rise,  even  when  the  Bz2  magnetic  pressure  is  confined  to  a 
subset,  or  layer,  of  field  lines.  Also  note  that  the  thickness  of  the  heated  layer 
will  increase,  just  as  the  field  lines  in  model  4  (Figure  2)  were  found  to  spread 
preferentially  in  the  layer  of  shear.  In  both  cases  the  effect  of  the  expansion  is 
to  alleviate  the  pressure  increase  brought  about  by  the  imposed  changes:  in 
the  magnetic  problem,  a  decrease  in  Bx2  +  By2  tends  to  offset  the  pressure 
increase  due  to  Bz2,  while  in  the  hydrostatic  problem,  a  decrease  in  particle 
density  offsets  the  pressure  increase  due  to  the  increase  in  temperature.  An 
important  point  is  that  the  expansion  is  entirely  upward.  There  is  no  external 
boundary  to  inhibit  this  upward  expansion,  whereas  a  downward  expansion  is 
prevented  by  the  rigid  photosphere  and  the  cap  on  the  bottom  of  the  tube. 

In  hydrostatic  equilibrium,  thermal  pressure  gradients  are  balanced  by 
gravitional  forces,  while  in  a  force-free  magnetic  field,  magnetic  pressure 
gradients  are  balanced  by  magnetic  tension  forces.  Analogous  to  equation  (5.1) 
therefore  is  the  result  that  the  magnetic  pressure  at  any  point  along  the  y-axis 
of  a  force-free  arcade  is  equal  to  the  integrated  effect  of  the  magnetic  tension  in 
the  overlying  field: 


— B(y)2=- 
871 


B, 


5B 

~cbTdy  ' 


(5.2) 


If  all  the  field  lines  of  an  arcade  rise,  then  the  field  strength  and  magnetic 
pressure  at  the  arcade  axis  (y  =  0)  must  certainly  decrease  (for  Bz  =  0  at  the 
axis;,  and  so  too  must  the  total  magnetic  tension  force  integrated  along  the  y- 
axis  from  zero  to  infinity. 
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In  Figure  7b  we  show  the  integrated  magnetic  tension  forces  in  the 
unsheared  arcade  of  Figure  2a  and  the  sheared  arcade  of  Figure  2b,  integrated 
along  the  y-axis  between  the  limits  of  0  and  y,  and  plotted  as  a  function  of  y. 
The  two  curves  indicate  the  cumulative  effect  of  magnetic  tension  as  one  move 
upward  along  the  y-axis  from  the  photosphere.  The  curves  clearly  diverge, 
with  the  curve  for  the  sheared  (inflated)  field  falling  below  that  of  the 
unsheared  field,  as  expected  from  equation  (5.2).  At  the  upper  boundary  of  the 
model  (y  =  60),  the  tension  integrals  differ  by  17.04%.  This  agrees  very  closely 
with  the  17.22%  difference  in  magnetic  pressures  at  the  arcade  axis, 
indicating  the  high  degree  of  accuracy  of  the  models  and  the  unimportance  of 
the  far-removed  upper  boundary.  In  models  where  the  upper  boundary  is 
important,  the  nonzero  pressure  at  that  boundary  must  be  added  to  the  right- 
hand-side  of  equation  (5.2)  (or  equation  [5.1]).  It  is  possible  for  the  inner  part  of 
a  sheared  field  to  become  depressed  in  such  a  model,  although  this  does  not 
seem  to  be  the  problem  with  the  Low  and  Nakagawa  result. 

By  comparing  Figures  7a  and  7b,  we  see  that  the  reduction  in  magnetic 
tension  in  the  sheared  arcade  occurs  mostly  in  the  region  of  greatest  shear, 
where  Bz  is  large.  We  can  express  the  tension  force  as 


F 


tension  ' 


471 


B 

R 


(5.3) 


where  R  is  the  local  radius  of  curvature,  given  by 


R(x=0)  =  - 


B*(dBy 

“bT"!  nr  I 


(5.4) 
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along  the  y-axis.  Since  the  tension  force  is  reduced  in  the  region  of  shear,  we 
know  that  the  magnetic  pressure  gradient  must  also  be  reduced  in  order  to 
have  force  balance.  This  suggests  that  B2  will  be  relatively  large  in  the  region 
of  shear,  or  at  least  comparable  to  its  potential  field  value,  which  we  have 
confirmed  to  be  the  case.  The  significant  reduction  in  magnetic  tension  must 
therefore  be  due  to  an  increase  in  the  radius  of  curvature. 

In  Figure  7c  we  have  plotted  the  radius  of  curvature  as  a  function  of 
height  along  the  y-axis.  We  see  that,  indeed,  the  radius  of  curvature  is 
increased  in  the  region  of  shear,  by  up  to  a  factor  of  two.  This  "flattening"  of 
the  field  lines  must  be  due  entirely  to  their  being  stretched  in  the  z-direction, 
since  the  field  lines  also  become  elevated,  which  tends  to  decrease  their  radius 
of  curvature.  The  latter  effect  is  especially  evident  for  the  lower  lying,  weakly 
sheared  field  lines,  for  which  the  radius  of  curvature  actually  decreases.  We 
note  that  the  radius  of  curvature  of  the  outermost  field  lines  is  probably 
affected  by  the  upper  boundary  and  should  therefore  be  treated  with  caution. 

VI.  CONCLUSION 

In  summary,  we  have  conjectured  that  any  footpoint  shear  imposed 
upon  a  simple  arcade  magnetic  field  will  cause  the  entire  field  to  become 
inflated--aZZ  of  the  field  lines  will  rise.  We  have  shown  this  to  be  true  for  a 
number  of  different  shear  profiles,  including  two  which  maximize  Bz  a  finite 
distance  away  from  the  arcade  axis.  Field  lines  tend  to  spread  apart  in  the 
region  of  enhanced  Bz  (as  also  found  by  Amari  et  al.  1989),  but  all  of  the  field 
lines  risc--the  higher  field  lines  simply  rise  more  than  the  lower  field  lines. 
We  have  also  demonstrated  that  the  original  depressed  field  result  of  Low  and 
Nakagawa  (1975)  is  incorrect,  and  that  the  sheared  field  actually  becomes 
inflated. 
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We  note  that  field  lines  can  become  temporarily  depressed  in  response  to 
a  sudden  or  rapidly  changing  shear  (rapid  compared  to  the  characteristic 
Alfven  travel  time)  (e.g.,  Wu  et  al.  1983),  but  this  is  a  property  of  configurations 
that  are  not  in  equilibrium,  and  we  are  here  concerned  with  configurations 
that  are  evolving  quasi-statically  and  at  all  times  well-approximated  by  a  static 
equilibrium. 

Tables  2  and  3  show  that  the  degree  to  which  fields  become  inflated  can 
be  substantial,  especially  in  the  outer  portions  of  the  configuration.  This 
suggests  that  the  shearing  of  magnetic  fields  may  be  an  important  process 
leading  to  the  opening  of  magnetic  field  lines.  As  closed  field  lines  are  brought 
to  greater  heights  in  the  solar  atmosphere,  they  are  more  likely  to  be  opened  up 
by  the  thermal  pressure  and  dynamic  forces  of  the  plasma  (e.g.,  Pneuman  and 
Kopp  1971).  Thus,  we  might  expect  shearing  motions  on  the  surface  of  the  Sun 
to  be  associated  with  an  increase  of  the  interplanetary  magnetic  flux  and 
possibly  an  increase  in  coronal  hole  area. 

Finally,  we  have  shown  that,  although  field  lines  rise,  they  can 
nonetheless  by  flatted  by  the  stretching  due  to  the  shearing  displacement.  This 
has  important  consequences  for  the  formation  of  solar  prominences,  since 
prominences  cannot  form  along  field  lines  that  are  well-rounded  (i.e.,  with  a 
small  radius  of  curvature)  and  concave  downward  (see  the  discussions  in 
Priest,  Hood,  and  Anzer  1989  and  Amari  et  al.  1989).  We  find  no  evidence  for  a 
dip  in  the  field-a  region  of  upward  concavity-in  agreement  with  the  recent 
result  of  Amari  et  al.  (1989). 

We  note  that  low-lying,  relatively  flat  field  lines  occur  naturally  above 
"weak  field  corridors"  of  sharply  reduced  flux  (see  Figures  2  and  7c).  in 
agreement  with  observation  (Klimchuk  1987,  1988).  It  is  in  these  corridors  that 
active  region  filaments  are  found  to  occur  (e.g.,  Klimchuk  1986,  1987).  The 
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degree  of  flatness  of  the  field  lines  can  be  increased  by  shearing  their 
footpoints,  near  the  edge  of  the  corridor,  and  indeed  corridor  fields  and 
filaments  are  often  observed  to  have  a  high  degree  of  shear  (e.g.,  Klimchuk 
1987;  Zirin  1972).  It  is  not  surprising,  therefore,  that  active  region  filaments 
occur  preferentially  in  weak  field  corridors  where  the  field  is  highly  sheared. 
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TABLE  1 


FIELD  LINE  ELEVATION  RELATIVE  TO  THE  POTENTIAL  FIELD 


Model 

Ay 

(a=0.8) 

Ay 

(a=0.5) 

Ay 

(a=0.2) 

Actual 

u/u0 

Predicted 

U/U0 

1 

3.19 

6.26 

12.47 

1.96 

1.87 

(219  %) 

(171%) 

(107%) 

2 

1.18 

4.07 

8.60 

1.55 

1.51 

(107%) 

(111%) 

(74%) 

3 

0.45 

1.50 

4.49 

1.24 

1.26 

(41%) 

(41%) 

(38%) 

4 

0.11 

1.85 

4.05 

1.20 

1.22 

(10%) 

(50%) 

(35%) 

5 

0.05 

0.49 

2.46 

1.07 

1.13 

(5%) 

(13%) 

(21%) 

TABLE  2 

FIELD  LINE  ELEVATION  RELATIVE  TO  THE  POTENTIAL  FIELD: 


EQUATION  (3.6)  SHEAR  FUNCTION 

(MODEL  4) 

7 

m  a  x 

Ay 

(«=(). 8) 

Ay 

(a=0.5) 

Ay 

(a=0.2) 

Actual 

U/TJ0 

Predicted 

U/U0 

2.0 

0.15 

0.55 

1.26 

1.09 

1.08 

(14%) 

(15%) 

(11%) 

6.0 

1.18 

4.07 

8.60 

1.55 

1.51 

<  107% ) 

(111%) 

(74%) 

10.0 

2.84 

8.86 

16.64 

2.03 

1.99 

(258%) 

(241%) 

(142%) 
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FIGURE  CAPTIONS 

Fig.  1. -Distributions  of  the  normal  component  of  the  photospheric  magnetic 
field,  Bv(y=0),  given  by  equation  (3.3)  and  the  footpoint  z-displacement,  y(y=0), 
given  by  equation  (3.5),  plotted  as  a  function  of  distance  from  the  arcade  axis. 
The  distributions  are  antisymmetric  about  the  axis,  and  only  the  positive  x- 
halves  are  shown.  Plotted  points  indicate  the  spatial  resolution  of  the  model 
calculations.  These  boundary  conditions  correspond  to  model  4  in  Table  1. 

Fig.  2. --Contours  of  constant  a  for  models  having  photospheric  boundary 
conditions  given  by  equations  (3.4)  and  (3.5):  (a)  the  potential  field  case  Zmax  = 
0;  (b)  the  sheared  case  Zmax  =  6  (model  4).  These  plots  represent  end-on  views 
of  the  arcades.  Please  note  that  the  domain  of  the  calculations  (60x60)  is  much 
larger  than  shown. 

Fig.  3. --Shear  functions  showing  the  distribution  of  footpoint  z-displacements 
as  a  function  of  distance  from  the  arcade  axis.  The  functions  are 
antisymmetric  about  the  axis,  and  only  the  positive  x-halves  are  shown. 
Plotted  points  indicate  the  spatial  resolution  of  the  model  calculations. 
Number  labels  refer  to  model  numbers  in  Table  1. 

Fig.  4. -Adopted  shear  function  for  our  computation  of  the  Low  and  Nakagawa 
model  (Case  II).  Points  are  from  Figure  3  of  Low  and  Nakagawa  (1975). 

Fig.  5. -Contours  of  constant  a  for  our  computation  of  the  Low  and  Nakagawa 
model.  The  boundary  conditions  on  a  are  given  by  equation  (4.1)  and  the 
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boundary  conditions  on  y  are  given  by:  (a)  y(x,0)  =  0  (the  potential  field);  or  (b) 
equation  (4.3). 

Fig.  6.— Hydrostatic  analogy  to  the  sheared  force-free  magnetic  field  problem: 
(a)  semi-infinite  vertical  tube  filled  with  an  isothermal,  gravitationally 
stratified  gas;  (b)  resulting  equilibrium  after  the  gas  has  responded 
adiabatically  to  a  temperature  increase  of  the  particles  in  some  layer  (gravity 
is  assumed  to  decrease  with  height). 

Fig.  7. -Magnetic  field  quantities  plotted  as  a  function  of  position  along  the  y- 
axis  (height)  for  the  unsheared  arcade  of  Figure  2a  (dashed  curves)  and  the 
sheared  arcade  of  Figure  2b  (solid  curves):  (a)  Bz  component  of  the  field;  (b) 
cumulative  magnetic  tension  force  integrated  from  the  base  (y  =  0)  to  height  y; 
(c)  local  radius  of  curvature  of  the  field.  The  units  of  (a)  and  (b)  are  arbitrary. 
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